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Abstract. Smectic liquid crystals are materials formed by stacking deformable, 
fluid layers. Though smectics prefer to have flat, uniformly-spaced layers, boundary 
conditions can impose curvature on the layers. Since the layer spacing and curvature 
are intertwined, the problem of finding minimal configurations for the layers becomes 
highly nontrivial. We discuss various topological and geometrical aspects of these 
materials and present recent progress on finding some exact layer configurations. 
We also exhibit connections to the study of certain embedded minimal surfaces and 
briefly summarize some important open problems. 

Keywords: mean curvature flow, foliations, liquid crystals 



1. Liquid Crystals: Phases in Between Liquid and Solid 

First characterized at the end of the 19th century, liquid crystals have 
gone from chemical curiosities to important optical and electronic de- 
vices (de Gennes and Prost, 1993). In fact, this article relied on liquid 
crystal technology to be produced! As their odd name implies, liquid 
crystals have some properties of a crystal and some of a liquid. Whether 
nematic, cholesteric, or smectic, all these phases flow (like a fluid), at 
least in some directions, while transmitting torques and forces (like a 
crystal) in others. The phase we consider here is the smectic-A phase. A 
crystal is characterized as a periodic structure (the simplest of which are 
lattices) . A d-dimensional crystal is a periodic structure in d dimensions 
with d independent periodicities. However, it is possible to have a lower 
number of periods than than number of dimensions. Focussing on d = 3 
dimensions, a one-dimensional periodic structure is a smectic, a two- 
dimensional periodic structure is a columnar phase, a three-dimensional 
structure is a crystal, and a phase with no periodicity is a fluid. 

In this article we will consider smectic phases and discuss minima 
of the free energy 

F= \l v d * X " !) 2 + ^i(V-N) 2 } (1) 

where the density of material is p{x) oc cos [27r<3?(x)/a], a is the one- 
dimensional spacing of the layers and N = V3>/|V3>| is a unit vector 
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field, normal to the level sets of <E>. The constants B and K\ are the 
bulk and bend modulus, respectively, and the integral is over the sample 
volume V. In particular we are interested in minima with non-trivial 
topology on domains which are not simply-connected. Before discussing 
those solutions, we will review the origins of the functional in (1). 

A one-dimensional periodic structure in three dimensions can be 
represented as an L2 summable density wave in R 3 

00 

p(x) = ( 2 ) 

j=-oo 

where p* = p-j, so that p(x) real, and q is a vector pointing along 
the direction of the modulation with magnitude 2ir/a, where a is the 
spacing between the "layers". Each Fourier mode pj for j 7^ is an 
order parameter for the smectic: if pj = for all j 7^ 0, then the 
material is a fluid, while if some pj is nonvanishing, it is a smectic. 
Since a is the spacing between the layers, it is sufficient to keep track of 
just pi since the presence of the fundamental mode is enough to signal 
the order. Truncating (1) to the first non-trivial mode and writing 
Pi = \pi\e~ 2mu / a (Chaikin and Lubensky, 1995) we have 

p(x) = po + 2\p 1 \ cos 

To study deformations and fluctuations of this periodic structure, we 
promote u from a constant to a scalar field u(x). This results in an 
Eulerian description of the material, i.e. the layers sit at surfaces of 
constant 3>(x) = q-x — u(x) or, equivalently, at the peaks of the density 
wave. An alternative description is Lagrangian, in which we label each 
fiducial plane defined by q ■ x n = n and consider deformations about 
these. When the deformations are small they can be represented as a 
height function or graph over each plane. If we choose our coordinates 
so that q = z, then 

n = h n (x,y) - u[x,y,h n (x,y)] (4) 

However, a unique solution for h n requires that the deformed surface 
not have any overhangs and, moreover, awkwardly mixes spatial co- 
ordinates with the Lagrangian coordinate n which labels the original 
surface. Despite this, as we will see subsequently, in restricted cases the 
Lagrangian representation is a useful tool. 

There are three essential actors in the free energy of smectics: the 
layer compression or strain, the bending or mean curvature, and the 
topology of the surfaces characterized by the intrinsic or Gaussian 



2tt 

— {q-x-u) 



(3) 
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curvature. Using the phase field <£, related to the density p via p oc 
cos (27r<I>/a), the unit layer normal is 



N 



(5) 



At a point on a surface, the two principle radii of curvature R\ and i?2 
determine the two curvatures Ki = \jR%. These two curvatures can be 
combined into the mean curvature H = \ (ki + K2) and the Gaussian 
curvature K = k\K2- The former is a measure of the bending of the 
layers in space, while the latter is an intrinsic quantity which is a local 
measure of geometric distortion and can be applied globally to gather 
topological information. Fortunately, these quantities may be expressed 
in terms of the layer normal: (Kamien, 2002) 



H = iv-N 



Q(x,y,z)=c 



K 



-V • [N (V • N) - (N • V) N] 



$(x,y,z)=c 



(6) 
(7) 



where the expressions are evaluated on one of the level sets of <F By 
working in this set of coordinates, it is possible to free these geometric 
quantities from the surfaces and, as we will show, allows us to relate 
different terms in the free energy. 

Because smectics have a preferred layer spacing, a, there is nec- 
essarily a strain energy associated with changing this spacing. While 
the two curvatures are natural from the geometric point of view, an 
expression for the compression strain is not as well defined. Returning 
to the family of equally spaced layers, we have na = &(x n + cmN) and 
so, differentiating with respect to n, find that 



N • V$(x n + anN) 



(8) 



which implies that |V$| = 1. Thus the strain, when expressed in 
terms of must vanish when this condition holds (Kleman, 1983). 
For instance, we can define the strain 



= 1 - |v$| = 1 - N • V$ 



(9) 



or an alternate strain u zz as 



^Z2 ~ 



1 - (V$)' 



(10) 



both of which vanish when |V<£| = 1. Moreover, in terms of u, the 
Eulerian displacement field, it follows that u zz « d z u + O (u 2 ) and 
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d z u + O (u 2 ) are identical to linear order in u and differ only 
in their nonlinear terms. Therefore, when considering the quadratic 
compression energy: 

F c = l -J d\Bu 2 zz (11) 

the harmonic term is independent of the specific choice of u zz - the 
linear term will always be d z u. While u zz arises naturally from a Landau 
theory description of the nematic to smectic transition (de Gennes, 
1972), there is no compelling reason to choose it over u zz - both are 
invariant under rotations and vanish for equally spaced layers. The con- 
stant B is known as the bulk modulus and controls compression modes 
of the layers. As we will show in the following, the use of u zz allows 
us to further develop the geometric theory of smectics 1 . In a similar 
vein, one could argue that the mean curvature H could be replaced by 
Iy 2 $. However, the nonlinear form H = • N unambiguously has 
the Laplacian as the first term in an expansion in powers of u, while 
there are many nonlinear generalizations of the Laplacian, including 
the full mean curvature. 

Putting this all together we have (Santangelo and Kamien, 2005a) 

F = — / d 3 x {[1 - N • V$] 2 + A 2 (V • N) 2 } , (12) 

where the last term is the contribution of the layer curvature to the 
energy, and is quadratic in the radii of curvature. The constant A = 
y/Ki/B has units of length and sets the natural lengthscale for defor- 
mations. For the sake of consistency, one should add to this expression 
a term of the form 

F K = K I d 3 xV ■ [(V • N)N - (N • V) N] , (13) 
Jv 

coming from the Gaussian curvature. Since this term is a total deriva- 
tive, it can be reduced to an integral on the boundary and, therefore, 
cannot play a role in the minimization of the energy that determines 
the configuration of the layers. The usual ground state which minimizes 
(12) depends on an arbitrary unit vector n: = n • x, corresponding 
to equally spaced layers along the n direction. The free energy of this 
solution vanishes since N = n is constant and |V$| = |n| = 1. 

It will become clear why we have chosen this nonlinear form for 
u zz . In the following section, however, we will linearize this theory - a 
process which results in an unambiguous form for u zz . 

1 Indeed u zz = 1 — — 2u zz = u zz + O (u 2 zz ) and so energies based on u zz will 
differ from those with u zz only by anharmonic terms. 
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Only in very few cases can we satisfy the terms in the energy individ- 
ually. The difficulty of this problem has to do with the nontrivial way 
that the terms in equation (12) interact. It is instructive to consider a 
family of evenly spaced layers x n (si, S2) parameterized by (si, S2), with 
x„ = xo + naN. At any point (s±, S2) on xo, we have the two principal 
curvatures ki 5 o and K2,o- The radii of curvature of the nth layer are just 
Ri,n = Ri,o + na from which it follows that 

1 + naKifi 

From the definitions of H and K it follows that the curvatures of the 
nth layer are: 

K ^0 ri( .s 

1 + 2naH + n 2 a 2 K { ' 

H - H ° + UaK ° (16) 
l + 2naH + n 2 a 2 K [ ' 

Thus if the layers are evenly spaced, the only way for the bending to 
vanish for all layers (H n = 0) is for Kq (and thus K n ) to vanish (Di- 
Donna and Kamien, 2003). This interplay between even spacing, the 
Gaussian curvature, and the mean curvature is a recurring motif in this 
work and will be discussed further in the following. It is clear that in 
order to have configurations with non-vanishing curvature, the layer 
spacing must necessarily vary. 

1.1. Topological Defects 

The ground state or equilibrium configuration of the field u subject to 
boundary conditions is the goal of our analysis. While the boundary 
conditions at infinity can usually be taken to be u = or Vu = 0, 
we may introduce additional boundary conditions via topological de- 
fects (Chaikin and Lubensky, 1995). In particular, because $ is a phase 
field, it can wind by multiples of 2ir around defects. In the so-called 
"London limit", we simply remove a line from M 3 around which the 
phase winds. It is possible to instead allow the density p\ to vary in 
space: the solution to the full equation has p\ vanishing at the core 
of the defect rendering the energy finite. This extra complication is 
tangential to our discussion and we will use the London theory in the 
following. 

Much of the physical understanding of defects comes from the lin- 
earized form of the theory or, equivalently, the quadratic energy (Kle- 
man, 1983; de Gennes and Prost, 1993): 

B ' " (3 z uf + A 2 (Viu) 1, (17) 



l 3 x 
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The first term is the linearized compression strain and measures the 
stretching along the layer direction, taken to be z. We define = 
xd x +ydy, the gradient in the plane perpendicular to z, and see that the 
second term is the Laplacian in that plane, which we recognize as the 
linearized form of the mean curvature, H. It is important to note that 
the original theory was quadratic in the strain u zz and the curvature H. 
There is no reason, physical or otherwise, that this truncation should 
be valid - it is used primarily because of its solvability. In addition, 
the ambiguity arising from the freedom to choose different forms of the 
nonlinear free energy does not arise in the linear theory. In order for 
the density to be single- valued, we have 

Vu-dl = na (18) 

r 

for some n G Z whenever the contour V goes around the defect. The 
combination —na = b is called the Burgers vector (Chaikin and Luben- 
sky, 1995). This can be promoted to differential form by introducing 
a defect density m(x) which is a vector field pointing along the defect 
line (de Gennes and Prost, 1993). The defect line is along the curve 
r(s) 6 M 3 with local unit tangent t(s) = dr/ds and we write 

m(x) = na/ 5 3 [x - r(s)] t(s) (19) 
Jo 

It follows that V • m = for defect lines that either run off to infinity 
or are closed loops. From Stoke's theorem, we have the local condition: 

VxVu = m (20) 

To find the ground state configuration of the smectic in the presence of 
defects, we write u as the sum of a regular part a and a singular part 
r so that diu = di<j + T{ with V x r = m (de Gennes and Prost, 1993). 
This is most easily done in Fourier space where 

r(q) = =^P^ (21) 

where we have chosen V • r = 0, absorbing any of the longitudinal 
parts of r into a. Having enforced the boundary conditions via the 
introduction of the singular part of u, we insert the expression for Vu 
into the free energy and minimize over the remaining variable, a. This 
results in 



6a 



{[d + AVja-*(l-AM)^™} (22) 



Kamien_Santangelo.tex; 2/02/2008; 3:38; p. 6 



Smectic Liquid Crystals 



7 



so that 



and 



a(q) = 



Fi = 



q z (1 - \ 2 q\) {q x m y - q y m x ) 
d\ \z ■ (q x m)\ 2 



BX 2 



I 



(27T)3 9 2 + A 2g4 



(23) 



(24) 



Because of the anisotropy of the smectic, we can have defects which 
are parallel to the special, z axis and defects which are perpendicular to 
it. The former are called screw dislocations while the latter are known as 
edge dislocations (Kleman, 1983). Defects can, of course, have a mixed 
character when they are oblique to the z axis. For an edge dislocation 
pointing along the y axis, m = nay5{x)5{z) or m{q) = (2ir)nay5(q y ). 
In order to calculate the energy we will need a long length cutoff L y /2 
for q y and a short distance cutoff g max = 7r/£ for q z . We have 



edge 



BX 2 n 2 a 2 , 

-^^Jdq x dq Zql + x2qt 
BL y \[\n 2 a 2 fi™* 

J <?m ax 

BL y n 2 a 2 y/ir [X 



5(q y = 0) 



/(/max \ 
dq z —j= 
' Qm ax \f H 



V2 



(25) 



so that the energy per unit length (F/L y ) is finite (de Gennes and Prost, 
1993; Kleman, 1983). For a screw dislocation, m = naz5(x)5{y) and the 
free energy vanishes meaning that the smooth part of u can "screen" 
the strain from the boundary at the defect. The spatially dependent 
solutions for the edge and screw dislocations are 



Vliedge 




2tt (x 2 + y 2 ) 



(26) 
(27) 



respectively. 



1.2. The Edge Again 

It is instructive to notice that the result for edge dislocation is simply 
the gradient of 



u{y,z) = - Sg n(z) b -^=J_^dte- t2 l 2 



(28) 
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where we will now switch to the Burgers vector b = —na. For z > 0, 



we have (d z — Xd% 



while for z < 0, (d z + Xd%) u = 0. There is a 



discontinuity between the upper and lower half planes which precisely 
enforces the condition (18). To see this, we take the clockwise path 
in the yz plane, starting at (y,z) = (— oo,0 _ ), going to (+oo,0~), to 
(+oo,0 + ) and back to (— oo,0 + ). We have 



Vu-di 



[it(-oo, + ) - u(oo, + )] + [u(oo, + ) - u(oo, - )] 



+ [u(oo,0~) - u(-oo,0 )] 



2 2 



(29) 



Using this result, we can calculate the energy directly by noting that 

2 



B f ( r 1 2 

F l = — d 3 x \ d z u t AVin ± 2A^nV 2 L u 
2 Jm I L J 



(30) 



D 

2 Jm 



d z u =f AV^n 



±BX f d 3 x (v± • (d z uV ± u) - \d z (V ± uf 
Jm I 2 

We take M to be all of IR 3 . In the upper half plane, we pick the minus 
sign, while in the lower half plane the plus sign so that the first integral 
over M vanishes and the second becomes a surface integral over dM, 
the xy plane. We see then that the energy is 



BXL 




dx exp 




(31) 



where, again, we have introduced a short-distance z cutoff, By choos- 
ing £' = £/(647r 2 ) the two expressions for the energy can be made to 
agree. The important point is that the energy has the same dependence 
on the physical parameters, b, B and A. 

The insight of Bogomoln'yi, Prasad, and Sommerfield, was to turn 
this calculation around: by rewriting the free energy as a perfect square 
plus a total derivative we can bound the energy from below simply 
through the boundary conditions (Bogomol'nyi, 1976; Prasad and Som- 
merfield, 1975). Moreover, when there are solutions of the lower order 
differential equation in the perfect square satisfying the appropriate 
boundary conditions, the energy arises solely from the boundary condi- 
tions. Such states are called BPS states and their energy is sometimes 
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called "topological" because of its sole dependence on the boundary- 
conditions. 

In the next section we consider the rotationally invariant model with 
this decomposition in mind. 



2. Rotational Invariance: From Here to Nonlinearity 

The quadratic approximation employed in the previous section has the 
drawback that it does not preserve the invariance of the free energy 
under rigid rotations of the layered structure. This rotation invariance 
is preserved by the higher order terms that would appear in an explicit 
expansion of the energy in the deformation u(x,y,z). As it happens, 
these nonlinear terms are crucial for understanding a number of physi- 
cal phenomena not captured by the quadratic approximation, yet they 
render the energy too complex for a straightforward analytical solution. 

2.1. BPS Solutions 

In the full free energy, there are also exact minima that satisfy a BPS 
flow equation (Santangelo and Kamien, 2003; Santangelo and Kamien, 
2005a). We start with the nonlinear theory (12): 

F= |y> x | [1 _ N . V< | )] 2 +A 2 (V . N) 2j 

Following the discussion in the last section, we complete the square 
yielding: 

B f 

F=— d 3 x [(1 - |V$| - AV • N) 2 + 2A(1 - |V$|)V • n] (32) 

While the first cross term is a total derivative (V • N), the second is 
not. However, we write 

|V$|V-N = [N(V-N)] • V$ 

= [N (V • N) - (N • V) N] • V$ (33) 

where the additional term vanishes because N 2 = 1 and N||V$. Re- 
casting this in terms of the geometric invariants, we have 

\V$\H = -V • {$ [N (V • N) - (N • V) N]} - K<5> (34) 

Thus, we can convert the cross term in (32) into a total derivative plus 
a purely geometric quantity (Santangelo and Kamien, 2005a): 



B 
2~ 



(Fx 



(1_|V$| - AV • N) 2 + AXQK 
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+2AV • {N — $ [N (V • N) — (N ■ V) N]} 



(35) 



If we restrict ourselves to systems for which K = 0, we again have 
the BPS decomposition. Fortunately, edge dislocations introduce layer 
bending in only one direction, and hence a single edge dislocation 
could be a BPS configuration, as could be a collection of parallel edge 
defects. Whether or not such configurations exist depends on whether 
all boundary conditions of the problem can be fulfilled by the lower 
order "BPS" equation: 



|V$| = AV-N 



(36) 



We have studied this equation numerically by switching to Lagrangian 
coordinates for the layers. In terms of layers described by an embedding 
X n for layers labeled by n, this can be rewritten as the flow equation 



N • d n X n = 



1 - 2XH 



(37) 



where we have used d^/dn = a. If we start with a layer at n = with 
vanishing Gaussian curvature, this flow generates new surfaces with 
vanishing Gaussian curvature, enabling us to find an exact solution 
for smectics with multiple parallel edge dislocations and to determine 
the energy due to bending the surrounding layers. The nonlinear BPS 
equation precludes, on the surface, superposition. However, in terms of 
u = (z — <&)/\ and the new variables a = xj\f\z and (3 = z/X, we find 
an expansion around (3 = oo (Santangelo and Kamien, 2005a): 



= dpu 



1 



1 



a (1 - dpuf d a u + (1 - dpuy - 2d l p u (d a u) 

(38) 



P 2(1 -dpuf 

-4 (i - dpu) d a ud a df3u + 2(1- dpuf a'nj + o ^ 

Thus, for large f3, we may write our solutions as a power series around 
(3 = oo. The lowest order, non-trivial equation for u = f(a) is indepen- 
dent of (3 and satisfies 



-\ad a f+(d a ff\+d 2 J = Q (39) 

Through the Hopf-Cole transformation S = exp{//2} we have 

ad a S + 2d 2 a S = 0, (40) 

a linear equation. 
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Converting back to y and z coordinates, we find that u approaches 
a similarity solution to the equation 



d z S 



\d 2 y S 



0, 



(41) 



for z > and the same equation with opposite relative sign for z < 0, in 
analogy with the linear case. The solution is (Santangelo and Kamien, 
2003; Brener and Marchenko, 1999) 



u(y, z) = 2Asgn(z) In 



1 + 



a fe/(4A) 



-I 



//\/4AR 



"* 2 / 2 dt 



(42) 



The maximum strain lies along the parabolas defined by y 2 
Since 



4A|z|. 

is constant for all z it fol- 



on these parabolas, u(y/4X\z\, z 
lows that equation 42 never approaches the linear solution along these 
parabolas even when z — > oo. In this sense, the linear theory never 
approaches the nonlinear theory. (Brener and Marchenko, 1999; San- 
tangelo and Kamien, 2003) 



3. Nonzero Gaussian Curvature 

One difficulty of this approach is that it is no longer exact when K ^ 0. 
The form of the BPS flow equation suggests, however, that we write 
the Euler-Lagrange equation for the extrema of the energy using T = 
1 — |V3>| — AV • N. In this language, we find that (Santangelo and 
Kamien, 2005a) 

V- (nT-A^— j =2XK, (43) 

where P T = I — N <g) N is the projection operator into the tangent 
plane of the layers. When K = 0, it is clear that T = is a solution to 
equation 43 as expected. In addition, there will be T / solutions cor- 
responding to boundary conditions on u that differ from those allowed 
by the BPS construction. 

Equation 43 has the form of a conservation law for T, in which 
Gaussian curvature behaves as a source for T. In general, this pair of 
equations is quite difficult to solve, and a full analysis of equation 43 
has not yet been undertaken. Nevertheless, there is one case with K ^ 
that has been solved exactly in the nonlinear theory: the screw disloca- 
tion. In the case of a screw dislocation, m = nazS(x)8(y) and so a(q) 
(from (23)) vanishes and 

u S cvew = -^-tan -1 (-') (44) 



2tt \x 
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in the quadratic theory. However, taking <E> = z — u in the full theory, 
we find (r 2 = x 2 + y 2 ) 

= [-by/r, bx/r, 2-rrr] 
^JW + (2vrr) 2 

T = 1 - ^1 + 62/(2^)2 (45) 

Note that V • N = because the helicoid is a minimal surface. Evaluat- 
ing (43), we see that the helicoid also is an extremum of the nonlinear 
theory. Multiple screw dislocations are currently under investigation 
(Kamien and Lubensky, 1999; Kamien, 1999; Santangelo and Kamien, 
2005b). 



4. Summary and Open Problems 

The role of nonlinear elasticity in smectic liquid crystals is still a nascent 
field, and there are a number of open problems that are relatively 
unexplored. 

Understanding the arrangement of layers in a particular smectic 
material is challenging precisely because of the interplay and frustration 
between the geometry of the layers and their spacing. Though this 
problem can be solved exactly in some special cases, it is difficult to 
construct even reasonable conjectures for the shape of the layers in 
more general situations. Nonetheless, the few known exact solutions 
available compel us to conclude that the harmonic approximation for 
the energy of smectic liquid crystals fails dramatically. 

For the special case of K = 0, exact solutions can be found that sat- 
isfy an equation of lower order than the Euler-Lagrange equation (San- 
tangelo and Kamien, 2003; Santangelo and Kamien, 2005a). These low 
order equations allow the determination of the shape of all the layers 
once the shape of a single layer has been specified. Furthermore, the 
energy of these configurations is reduced to a topological term, which 
can be evaluated on the layers near the defect core. 

Some analytical progress for K ^ surfaces has also been made by 
directly summing the phase fields of screw dislocations. In the case of 
a twist-grain boundary, where the screw dislocations lie along the line 
x = 0, there is an unexpected connection to minimal surfaces (Kamien 
and Lubensky, 1999; Santangelo and Kamien, 2005b). Neither sums 
of screw dislocations, nor their minimal surface counterparts, however, 
are exact minima of the smectic energy - the true minimum likely lies 
somewhere in between (Duque and Schick, 2000). More generally we 
can ask if there are solutions with K ^ with energy determined 
entirely by the boundary conditions? 



Kamien_Santangelo.tex; 2/02/2008; 3:38; p. 12 



Smectic Liquid Crystals 



13 



An analysis of the BPS flow for T in the presence of Gaussian cur- 
vature has not yet been undertaken. How do the layers evolve in the 
presence of Gaussian curvature? Even when K = 0, there are nontrivial 
aspects of the layer structure. When there are two edge dislocations 
separated by several intervening layers, what is the structure of the 
layers in between and how do those edge dislocations interact? In other 
words, can the BPS decomposition shown here be extended to solutions 
with constant or varying Kl 

In addition, because both screw and edge dislocations are topolog- 
ically identical, there exists a family of dislocations with the same 
topology which can be found by rotating the defect core from being 
parallel (edge) to perpendicular (screw) to the layers. These defects will 
have both edge and screw components, and we expect that nonlinear 
elasticity will be important in describing both their energetics (as is the 
case in screw dislocations) and the geometry of the smectic layers (as 
in edge dislocations). Thus far, these defects have resisted theoretical 
descriptions beyond the harmonic approximation. What is the structure 
and energetics of defects with both screw and edge components? 

To summarize some open questions: 

• At long distances the extremals of (1) satisfied (39) which is actu- 
ally a linear equation. Thus, in this limit, we could find solutions 
for an arbitrary number of parallel edge dislocations at the same 
z coordinate by setting the discontinuity between two solutions 
of the diffusion equation, one diffusing "up" the other "down" 
(Santangelo and Kamien, 2003). ^,Can the BPS decomposition be 
used to find solutions for parallel edge dislocations separated along 
the z direction? 

• Single screw dislocations are helicoids and thus H = 0, suggesting 
solutions with vanishing compression, i.e. u zz = 0. However, u zz / 
for single screw dislocations and so it is unlikely that there are 
superposition solutions which satisfy u zz = 0. We are free to add 



with A 2 = K/B. lis there a BPS decomposition of this functional 
for H = 0? lis there a decomposition if we use a different power 
of K in (46), e.g. K 2 ? 

• <iAre there extremals of (1) which are not parallel to the z axis 
(screw dislocations) or lying in the xy-plane (edge dislocations)? 
^,How does the BPS edge solution deform into the minimal surface 
screw dislocation? 



F K (13) to (1): 




(46) 
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The solution of these problems has the added benefit that they can 
be verified directly through experiments on liquid crystals (Ishikawa 
and Lavrentovich, 1999) and the physical validity of any approxima- 
tions can thus be assessed. Despite the long history of theory and 
experiment on these layered systems, the nonlinear aspects of these 
problems are only now being explored. 
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